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Abstract 

The degenerate free Fermi gas coupled to a random potential is used 
to compute a.c. conductivity in various dimensions. We first formally di- 
agonalise the hamiltonian using an appropriate basis that is a functional 
of the disorder potential. Then we compute the a.c. conductivity at zero 
temperature using the Kubo formula. This a.c. conductivity is a func- 
tional of the disordered potential. The wavefunction of extended states is 
written as exponential of the logarithm. We use the cumulant expansion 
to compute the disordered averaged a.c. conductivity for Gaussian disor- 
der. The formula is valid if a certain linearization approximation is valid 
in the long- wavelength limit. 

1 Introduction 

A simple theory of the Anderson transition is presented that directly computes 
measurable quantities such as a.c. conductivity. The relevant literature on 
this subject is vast and we shall not attempt to be exhaustive in surveying 
it. Anderson's pioneering work on localization[l] was followed by the work of 
Abrahams et.al. [2] and later on a more rigorous formulation of the notion of 
disorder averaging was given by McKane and Stone [3]. This relates to a single 
electron in a disordered potential. The classic review of Lee and Ramakrishnan 
[4] includes many references on the literature concerning the degenerate electron 
gas in a disordered potential. A more recent review is by Abrahams and [5]. 
Belitz and Kirkpatrick[6] have a review that is aimed at theorists. The present 
write-up is intended to give a simple derivation of the Anderson transition in 
three dimensions. 

The basic idea of this work is as follows. We formally diagonalise the hamil- 
tonian of a single electron in a fixed disorder potential. The wavefuntion of the 
electron then may be formally written as, 

^(x)^e*W^(x) (1) 
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Here d is an appropriate normalization constant. The novel feature involves 
rewriting the square of the amplitude i?i(x) also as an exponential. 

ifc(x) - e A * (x) (2) 

Then one uses the observation that only extented states participate in the con- 
duction. Thus we may write, 

&(x) = ^ e ,k -^ i(x) e5 A ' W (3) 
vV 

Then we make the ansatz that < 9i >di S =< Aj >dis~ 0. This ansatz is shown 
to be consistent in the main text. Then one uses the Kubo formula to eval- 
uate the a.c. conductivity. One simply plugs in the form in Eq.( 3) into this 
formula and uses the cumulant expansion to evluate the disorder averaged con- 
ductivity. In the cumulant expansion, we encounter correlation functions such 
as < 6i(x.)6j(x ) >dis, < #i(x)Aj(x ) >dis and so on. These have been explicitly 
computed in the appendix. It can be shown that in some approximate sense, the 
higher order correlation functions all vanish. Hence we may deduce a formula 
for the disorder averaged a.c. conductivity provided it is legitimate to use only 
the linearization approximation. Our formalism does not have an arbitrarily 
chosen cutoff. A natural smooth cutoff emerges from a careful treatment of cur- 
vature effects by retaining double derivatives, in other words, by not linearizing 
the bare fermion dispersion. In passing we note that this formalism involves no 
bosonization or other advanced field theoretic ideas. It only involves the simple 
Kubo formula. The answer for the a.c. conductivity may be written down as 
follows for Gaussian disorder. It has been reduced to quadratures. 
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(4) 



where, 



(JjiiUdis) ■ JaiUdi.)) =-\\ ddx e^- k ^ x e5 <£ ' 2> F^k,) (5) 



1 V2 ^_ m 2 A 2 v / 1 1 \ 2 m 2 A 2 v / 1_ , 1 ^ 2 



<E Z >- 



(6) 

F(k- 1c) - 1 V 2m2A2q2 c -fa- 
1 ^ 2m 2 A 2 ... 1 v^/ 2m 2 A 2 .. . 2m 2 A 2 .. \ _ 4Q x . 

1 ^ 2m 2 A 2 , . N 1 ^/ 2m 2 A 2 , . ^ 2m 2 A 2 . \ _ jq X1 
' ^ + 2V E A^q) (- q ) + 2F E { A (-q) ^(q) ^ + 75?^ J 6 " 1 

(7) 

A(q) = -k i .q+ q r (8) 

The presence of the q 2 /2 makes the integrals hnite at high momenta. This comes 
about by retaining double derivaties, in other words by taking into account the 
parabolic nature of the band. The remaining sections explain how these formulas 
are derived. Before we do this we would like to point out some general facts. 
Since, 

/(ki.kj-) = (JjiiUdis) ■ 3ij{U dia )) = (\J t3 (U dls )\ 2 ) > (9) 
We may sum over the angles first to write, 

1/(11^1,1^1) = J2f(^o)>0 (10) 
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/ 2\ i i-co r-co £,2 fc2 ,2 £.2 

(11) 

For tyllu <C hp, we have, 

ite[tr(w)]~ />f(1 + fc F ) (12) 

One really important question that we have to face is whether the zero fre- 
quency limit of the a.c. conductivity (obtained via Kubo formula) is the d.c. 
conductivity. According to the review by Lee and Ramakrishnan[4] this is true 
only in three dimensions. In one and two dimensions, it appears that according 
to the results presented, the d.c. conductivity may not be obtained by tak- 
ing the zero frequency limit of the a.c. conductivity. They show that the a.c. 
conductivity diverges logarithimically for small frequencies in two dimensions 
whereas the d.c. conductivity is strictly zero. A frequency cutoff emerges in 
these dimensions that tells us that for frequencies smaller than this cutoff the 
formulas given there break down. According to the analysis given in Appendix 
C of this article, there should be no problem in any number of dimensions. The 
zero frequency limit of the a.c. conductivity is in fact the d.c. conductivity. In 
three dimensions we may expect to find the mobility edge while computing the 
d.c. conductivity. 

^ d.c. ~ f( k F, kp) (13) 

Results So Far : 

So far the results have been very disappointing. In one dimension, we expect 
to see that the d.c. conductivity is zero. Instead we find that it is not although 
there appear to be many subtle cancellations. It is likely that the linearization 
is not adequate and retaining some nonlinear terms may be required. The main 
culprit seems to be an approximation that replaces a certain variable quadratic 
in the 9 and A by it disorder average in an effort to render the equations linear 
differential equations. However when evaluated this average is zero and hence it 
is a really bad idea. Hwoever the author had hoped that it would not be serious. 
But it is. Future attempts may try and address these problems, although there 
are many good ideas in this preprint. 

2 Some Technical Musings 

It appears that the mathematical literature on the subject of quantum par- 
ticles in random potentials is vast [7]. It is possible, indeed likely that many 
mathematically rigorous results are known regarding this problem. But this 
does not prevent the authors from making some remarks that more knowledge- 
able readers may choose to critique. In particular, the author is uncomfortable 
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with the notion of disorder averaging. Nature chooses its potentials based on 
the distribution of impurities, defects and so on. This potential is fixed and 
well-defined for a particular distribution of these imperfections. The physicists' 
ignorance of the precise nature of this potential is not a license to average over 
these potentials. Nature does not average, people do. But are people justified in 
averaging ? In other words can averaging simplify the problem without washing 
out essential physics ? In order to answer this question we have to make the 
following conjectures. 

DefnO : Let Ltd be the set of all potentials U(x) in a fixed spatial dimension d. 
Defnl : Let Td be the set of all potentials U (x) in a fixed spatial dimension d 
that has the following property. They all lead to the same exponent S for the 
frequency dependence of the a.c. conductivity. In other words, each of these 
potentials predicits that Re[a(co)] <~ u s (in some region of ui with possibly some 
additive part independent of lo) with the same 5. 
Conjecture 1 : Td is dense in hid- 

If Conjecturel is valid, then one may average over all these 'sufficiently er- 
ratic' potentials and expect to extract <5 which is all that physicists care about. It 
is possible that 5 may be extracted from a numerical solution of the Schrodinger 
equation using a specific U that belongs to the set Td- But this would involve 
using the computer for more than checking one's email, and not everyone likes 
that. 

Defn2 : Let be the set of all potentials C/(x) in spatial dimension d = 3 
that has the following property They all lead to the same exponent (3 for the 
mobility edge exponent. In other words, each of these potentials predicit that 
Od.c. ~ (Ef — Ec) 13 9(Ep — E c ) with the same (3. However for different potentials, 
E c - the mobility edge, may be different. 
Conjecture 2 : JF 3 is dense in U 3 . 

If Conjecture2 is valid, then one may average over all these 'sufficiently 
erratic' potentials and expect to extract (3. 

Thus the validity of the process of averaging over potentials rests crucially 
it seems, on all these sufficiently erratic potentials predicting the same expo- 
nents and on these sufficiently erratic potentials spanning nearly all possible 
potentials. 

If both these are satistifed then one may average over all potentials and 
extract the exponents, or, if one is better at programming, choose a particular 
potential from this set, numerically solve the Schrodinger equation and extract 
the exponent from there. In either case we should get the same answer. A final 
conjecture seems appropriate. 

Conjecture 3 : Let M 3 have an exponent (3 and T d have an exponent 5 , 
then (3 = (3 and 5 — 8 . In other words, these exponents are unique. 

With powerful computers now available, purely analytical methods such as 
this work may seem passe, but a closed formula for the a.c. conductivity that 
one can stare at (and one that is hopefully right) and admire has a charm that 
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a cold data file on the hard disk is unable to duplicate. Besides, with Coulomb 
interaction, the problem becomes intractable numerically, however, one may 
expect to combine the sea-boson method with the present one to extract the 
exponents analytically. 



3 A.C. Conductivity Using Kubo Formula 

In this section, we derive a formula for the a.c. conductivity in terms of the 
total momentum-momentum correlation function. To derive this, observe that 
the relevant hamiltonian that couples to external fields is of the form, 

H ext (t) = -M A ext (t)-P (14) 
m 

where P = k c^Ck is the total momentum operator. Thus we may define 
formally the net momentum in the presence of the external field in the imaginary 
time formalism in the interaction representation as 



T S P(t) 

(p W) . (15) 

where the S-matrix is defined as, 

S = e l V '*£ A «*(*)-P(t) (i 6 ) 

Here P(i) evolves according to the time-independent hamiltonian. The net 
current density in the system is given by J(t) = (|e|/my) P(t), this is in units 
of charge flowing per unit area per unit time in three space dimensions. From 
Ohm's law we expect, 

r if> > 

(J(t)) - / dt a(t-t)E ext (t) (17) 
Jo 

The a.c. conductivity for complex frequencies is then given by, 

-i/3 



a(iz) = 

Jo 



dt a(t) e~ z * (18) 



The d.c. conductivity is then given by, 

<Td.c. = + e) (19) 

where, e = + . 

»(«-*') (») 

E ext =0 



mV \5E ext (t') 
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* {t - t)= riV SE^t) (TS) (21) 



Since E ext (f) = -dA ext {t)/dt, we have, 



6Aext ^l=-e(t-t) (22) 



SEexti?) 

s „ J . r*' 



-S = S\-i[ dti 6{t -ti)^P(ti) ] (23) 
) V 7 m 



<*Eext(*') 

2 /— */3 



a(t-t') = - 2f 



m 2y 



jf <ftifl(t'-ti)^T<JP(ti)-<JP(t)^ (24) 



T£P(ii) ■ £P(t)^ =6(t!-t) 5^(k-k')JV(k,ti-t;k',0)+e(t-ti) 5^(k-k')JV(k,t-ti;k',0) 

kk' kk' 

(25) 

If we write for the dynamical number-number correlation function, 
iV(k,i;k',0) = (n k (t)n k /(0))-<n k (t)) (n k -(0)) = E e'^^'A^k, Ci , e j; k', 0) 

(26) 

where rik = c^Ck- We may write, 

(rapfa) • sp(t)) = e(ti - t) E( k ■ k') E 7v(k, c<> e j; k', oe-^-^)^-*) 

kk' ij 

+9{t-h) ^(k-k')^iV(k, ei , e3 ;k',0)e-^^)( t -* 1 ) (27) 

kk' *,3 

From the above formulas retaining only the terms that do not violate causality 
and using iz n — > u) — «0 + , we have at absolute zero, 

Re[a(u,U dis )\ = (^-) -^^(k.k) N(k, e h e r , k', 0)5 {u - a + ej) (28) 

kk' J J 

To compute this we have to first evaluate the dynamical number-number corre- 
lation function. Consider the hamiltonian, 

h = E + E ^# E 4 +q/2Ck - q/2 (29) 

k q V V k 

This may be diagonalised by the following formal transformation. 

Ck = E^( k ) di ( 30 ) 
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<c] c c k >=^|^(k)| 2 0(e F -e J ) (31) 

i 

n k (t) - J2 VjW^dfce-*'-** (32) 
y 

n k '(0) = E^'( k ')Vi'( k ')4 d *' (33) 

i' 3 

N(k, t; k',0) = ^(k)^(k)^(k')^(k>(e F - - e^e"^"^' (34) 

This means, 

2(k.k') JV(k,e i ,e J -;k',0) - J^(C/ dls ) • Jy^iaMe* - e,)0(e 4 - e F ) (35) 

kk' 

Here ^ (x) is the solution to the equation below, and <pi (k) is its Fourier trans- 
form. 

V 2 



+ U du (x)j 0i(x) - e, 4(x) (36) 
J«(tf«K.) - E k # ( k M"( k ) = - f d d x #(x)*V x ^-(x) (37) 



k 

i?e[cr(w, t/ dis )] = J ^ X! J ji( U <iis)-Jij{Udis) 9{e F -e j )e{uj+e j -e F ) 5{w-ei+ej) 

(38) 

Till now the discussion has been at the formal level. Now we would like to 
compute the disorder averaged conductivity assuming a Gaussian disorder. In 
other words, the averages over the disordered potential have to be performed 
using the conditions, 

< U dis ( X ) >= (39) 
< U dis {x)U dis {*) >= A 2 S d (x - x ) (40) 

Thus the sum over all configuarations of disordered potential keeps the sum 
/ d d x U dis (x) fixed and sums the conductivity obtained from each configuration 
and then computes the average. Several observations may be made regarding 
this. First, the energy e, for each choice of U d i S may be discrete and negative 
(bound state) or positive and continuous, corresponding to Bloch waves. The 
delta- function forces u) to be equal to the difference — ej. For small u), the 
difference e, — ej is likely to be comparable to u> only for the Bloch states. That 
is to say, that only the extended states participate in conduction. Thus we 
may polar decompose the wavefunction fa and assume that the magnitude is 
(roughly) independent of position, a feature characteristic of extended states. 
We have found that these simplifying assumptions though natural lead to some 
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divergences in the ultraviolet regime. Thus we shall have to retain the double 
derivatives. This is done in the appendix and the correlations between the A's 
and the 9's(see below) and amongst themselves are derived. Thus we write, 



0i(x) =e w '< x >fl?(x) 



(41) 



For some phase #j and Ri that is a function of the disordered potential. These 
obey the following set of coupled equations. 



V 2 & V0; • Vi?« 



2m 



2mR l 



= 



(V0j 



2m 



1 

8m 



2V 2 R t 
Ri 



(vgO 



Now we set, 



i2i(x) - ^ e A ^ x ) 



(42) 
(43) 

(44) 



Here C° is a normalization constant. It may be determined as follows. 



J d d x \Mx)f) = / d d x (Ri(x)) dis = Jd d 



Thus 



,- (e A,(x)^ =C O e |(A?(x)) = 1 
(45) 

(46) 



It will be shown in the appendices that {^A 2 (x)y is independent of position. 

Further since only extended states participate in the conduction, we may write 
#i(x) = ki -x+#j(x). In the appendices, an approximate scheme is written down 
that allows for a relatively simple computation of the average < 3{j (Udis) ■ Jji (Ud 
From Eq.( 38) it is clear that we would very much like to consider the eigenener- 
gies 6i as being independent oiUdis in the sense that we may replace ej =< e, >dis 
This requires some justification. Also from the appendix we see that 



,)> 



V 2 ^ + k t • VA, + W9 t • VA, = 



+ 



(W,) 2 + 



1 

8m 



2V 2 A 4 + (VA,) : 



= e» — C^di 



(47) 



(48) 



2m 2m 

The above Eq.( 47) and Eq.( 48) are absolutely exact. The approximations arise 
when we decide to linearize the above nonlinear partial differential equations 
by replacing the quadratic parts with their disorder averages. Obviously, this 
is justifiable only if we show that the fluctuations of the operators we have 
averaged out arc small, preferably zero. First we try and justify e, =< e» >dis- 
For this we have to evaluate < (e,— < ej >dis) 2 >^ and show that it is small 
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compared to < et >dis- To do this we take recourse to perturbation theory. 
Using elementary perturbation theory upto second order we may write, 



= S + / it ^ (x) + h E -kzk 1/ ddx ei{ki - ki) ~ u ^ 



K i 2m 2m 

Now, 



2 (49) 



k 2 A 2 v-^ 1 

< « >d ™= 2^ + -^ Tf~kf ^ 
k j/ k i 2m 2m 



Therefore, 



2 2 A 4 \ - 1 1 A 4 ^ 

< >di s - < ei > dls ~ -JT+TH 2^ 7k 5 k?7 7k? f2k,-k,)n + T?2 1^ 



2m 2m J \2m 2m y k j ^2m 

(51) 

The sums over kj in Eq.( 51) are finite in all three dimensions and hence the 
fluctuation e; is vanishingly small in the thermodynamic limit. Even other- 
wise, within this scheme the sum in Eq.( 50) diverges, hence the fluctuation 
(< e 2 >dis — < f-i >di S )^ is bound to be small compared to the mean. The 
proof that we may legitimately linearize the equations above (Eq.( 47) and 
Eq.( 48)) is given below. First we note that two random variables Oi(x) and 
2 (x) may be replaced by their means if, 

< O l (x)0 J (x') >=< 0;(x) >< Oj-(x') > (52) 
for i, j = 1,2. Thus we would like to make the following identifications. 



\2m 



V6i ■ VA 4 w ( V0i • VA 4 ) = {from appendix) (53) 

\ / dis 

l ( VA ( ,^i(^---L<(VA ( ,s> J „ (54) 
2m 8m 2m 8m 

The first condition is satified if we ensure, 

' (Vxflj • V x Ai) (v x * 0j • V x - A,) ) « (55) 
In other words, 

V™V^, < 0i(x)0,-(x') > dis V™V^, < A^A^x') > dis 

+V™V™, < ^(x)A,(x) > dis V™V£ < A^x^x') > dis « (56) 
The second condition is obeyed if we ensure, 

[25{Vj t f i<5(V x A,) 2 ] [2<5(V X ^,) 2 - ^(V^A,) 2 ] 
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= 8 (V™ V£ (^(x)^(x ))) (V™ V£ (^(x)^-(x ))) 

+ ^ ( V x K' (A t (x)A,(x ))) (v?V£ (A 4 (x)A,(x ))) 

-2 (V?V£ (A 4 (x)^(x ))) (v™V£ (A 4 (x)^(x ))) 

-2 (v?V£ (^(x)A,(x ))) (v™V£ (^(x)A,(x ))) « (57) 
Finally, we have to also ensure that the cross correlation functions are zero. 

i (V™V^ < ^(x)^(x') >) (V"V^ (A 4 (x)^(x ))) 

-4^ ( V ™ V£ < ^(x)A,(x ) >) (v™V£ (A,(x)A,(x ))) « (58) 

The three Eq.( 56), Eq.( 57) and Eq.( 58) represent correlation functions of 
random variables that we have assumed may be replaced by their averages. 
Unfortunately, we have found that these conditions are too severe to be obeyed. 
Thus we have to be content at saying that the formulas are not exact, but 
are based on a plausible linearization assumption. This assumption is 'self- 
consistent' in the sense that no obvious inconsistency shows up. A hand-waving 
justification for at least Eq.( 54) is that the right hand side diverges when 
evaluated thus the fluctuation of this quantity being finite is certainly small 
compared to the mean which is formally infinite. This infinity may be absorbed 
by a suitable redefinition of the Fermi energy. On the other hand, Eq.( 53) is 
impossible to justify except that making it renders an elegant analytical solution 
possible since the equations are now linear. 
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4 Appendix A 



= X7Ln[Ri\ = VAi 



We then have to solve, 



(V0 4 ) 2 



w 2 e t + V0< • va, = o 



1 



2m 8m 
We now decompose 9i as follows 



2V 2 A, + (VAO 2 



= e» - ?/<«. 



6i(x) = k, • x + ^(x) 
The reduced system may be written as follows. 

V 2 9 t + k t • VA, + W 4 • VA, = 



(59) 

(60) 
(61) 

(62) 
(63) 



k 2 | (W t ) 2 | kj^Wj 1 
2m 2m m 8m 



2V 2 A, + (VA l ) 2 = ei -C/ dis (64) 



Now we make the assertion that < 9i >= and < Aj >= 0. By taking the 
average of the above equations we may deduce that 



V9i • VA, ) = 



(65) 
(66) 

(67) 
(68) 

m 4m 

The above two equations may be used to iteratively compute the various corre- 
lation functions. First we have, 



2m + 2m 8m V VA,J ' ' 

From this we may also deduce, 

V 2 ^ + k, • VA, w 

k, • Vh V 2 A, 

-U di. 



V 2 < ^(x)C/ dis (x ) > +ki • V < A,(x)J7 djs (x ) 



(69) 



kj ■ V < 6» i (x)t/ d »(x ) > V 2 < A i (x)t/ d »(x ) > 
m 4m 

We solve this by Fourier transforms. 

1 



< U dis (x.)U dis (x) >= -A 2 <5 d (x-x ) 
(70) 



< ^(x)C/ d2S (x ) >= _£Vio(qi) e-^-( x - x ) 



(71) 
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< Ai(x)U dia (x) >= ^^F 20 (qi) e -*i-(*-*') 

Thus we have, 



V 
q 



-q 2 F 10 (qi) - * (ki ■ q)F 2 o(q«) = 

1 2 

_ z ^_H Flo(qj ) + -S-F 20 (qz) = - A 2 
V 2 < 6' J (x)A J (x') > +k 4 • V < A i (x)A j (x') > w 



k 4 • V < 6> i (x)A j (x ) > V 2 < A i (x)A j (x ) > 
m 4m 



< [/ dls (x)A i (x ) > 



< ^(x)A,(x ) >= I^F 12 (q;ij) e -M*-* ') 



q 



< A 4 (x)A,(x ) >=^F 23 (q;ij) e "W.(x-x') 



q 



Fi2(q;y)= (-^ + ^^) H(k,-q)F 20 (-q 7 )] 

F 22 (q;ii) = (-^ + 55^) 1 [q 2 F 20 (-qj)] 

V 2 < ^(x)^-(x') > +k 4 • V < Aj(x)0,-(x') > m 
k t .V<^-(x)> _ v 2 <A^(x)> m _ < § . , > 

m 4m 

< ^(x)^(x ) >=^fu(q;y) e -*q-(*-*') 



q 



< A t (x)^-(x ) >= I^F 21 (q;y) e -**(x-*') 
q 

Jii(q;y) = f-J- + ^^-) [-i(ki • q)F 10 (-qj)] 
V 4m m / 
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-1 



^2i(q;y) = (-^ + i^FM-qj)} (88) 

From the above equations it is clear that, 

< V0 t • VA, >= (V x • V x 0| x=x ' < 0i(x)Ai(x' ) > 

q 

e *-2^-w^ q (90) 

q x 7 

In one dimension, we have, 

kf mA 2 f°° / 1 1 \ 

e » = + tt^- / ; r- (91) 

2to 2Trki J q \x — Kj x + ki J 

If we interpret the above integral as the principal part then we have, 

k 2 



(92) 



In two space dimensions, it appears that we have to be more careful. In par- 
ticular, we have to introduce a large momentum cutoff that may not be easily 
dropped. 



k 2 mA 2 dx 



We take the point of view that this may be absorbed by a suitable redefinition 
of the Fermi energy. Thus in all three dimensions, we take the liberty to set 
€i = k 2 /2m. 
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5 Appendix B 

Thus, 

(JijiUdis) ■ JjiiUdis)) = - J d d x J d d y S d (y - x) J d d x J d d y S d (y' - x') 
(V y -V y O (#(x)&(y)#(x')&(y')) 

= J ddx j dd v sd (y - x ) J ddx ' J dd v sd (y' - x ) 

(V y -V ') e * k i-(y'- x ) e ik i-(y- x ') / e ^ 3 (y) + ^Ai(y) el 0i(y') + ^A I (y') e - 4 e,(x) + iA s (x) e - 4 6 J (x') + lA 3 (x') 



Thus we have, 

^■(y) + p J (y) e4 e;(y') + iA 1 (y') e -i0 1 (x) + lA 1 (x) e - 4 e j (x') + iA 3 (x )\ = e i<£ 2 > 



< ^ 2 >= ((^(y) + ^(y) + Wi(y') + ^(y) - i0» + ±A<(x) - ^(x) + ^(x')) 2 

= < ef(y) > < A 2 (y) > - < eUy) >+\< A 2 (y ) > 

- < fi?(x) > +1 < A 2 (x) > - < 2 (x') > +\ < A 2 (x ) > 
+^ < ^•(y)A j (y) > -2 < ^(y)^(y') > +i < ^(yjA^y) > +2 < ^(y)^(x) > +i < ^(y)A,(x) > 
+2 < ^(y)^(x') > +* < ^-(y)A,(x ) > +i < A,( y )^(y ) > +1 < A,(v)A,(y ) > 



-i < A,(y)^(x) > +1 < A J -(y)A i (x) > -i < A,(y)^-(x') >+\< A J (y)A,(x') > 
+^ < ^(y')Ai(y') > +2 < ^(y')^(x) > +i < i (y')A i (x) > +2 < ^(y')^(x') > +^ < ^(y')A,(x') > 
-i < Ai(y')fl>i(x) > +i < A i (y')A < (x) > -i < Aj(y')0,-(x') > +1 < A^y'jA^x') > 

-i < l (x)A l (x) > -2 < ^(x)0j(x') > -i < ^(x)Aj(x') > 

1 

2 



-i < Aj(x)0,(x ) > +- < Ai(x)Aj(x )>-i< 6j(x )A j (x ) > 



<£ 2 >= J B (i,i)+i^A 22 ,(q;i,i) e-^^V-^A^q^j) e^'-x) 



q q 



l£A 21 ,(q;i,j) e -^(y-') + I £ A, v (q;M) e^*'-*') 



q q 
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+^A 21 (^i,j) e-My-^ + JL^^q.^-) e -iq.(x-x') (g4) 
q q 

Eo(iJ) = -\ Yl F ii(q;ii)+^7 E ^(q; jj)-^ E F n(q; ^(q; ll ) 

q q q q 

(95) 

A 22 '(q;i,j) = -2 Fn(q;ji) +* Fi 2 (q;ji) +i f2i(q; j») + ^ F 2 2 (q; j'i) (96) 

4j'i(<i;m) = 2 ^n(q;") + * ^(q;**)-* *2i(q;i*) + ^ ^(q;**) (97) 
^21' (q; *>i) = 2 ^n(q; .?'.?) + * ^(q; ±0 - * F 21 {q-jj) + x - F 22 ( q; jj) (98) 

A^(q;i, j) - 2 F u {q;ij)+i F 12 (q;ij) - i F 21 (q:ij) + 1 F 22 (q;ij) (99) 
^2i(q;», j) - 2 fii(q;j'i) +i Fi 2 (q; ji) - i F 2 i(q;.7*) + \ F 22 (q;jt) (100) 

Ai' (q; *>i) = - 2 F n(q; y) - * *i2 (q; y) - * ^21 (q; y) + ^ ^22 (q; y) (101) 

(3 t3 (U dts ) ■ Julian)) - — ^r 1 - y ^ y dV e^- k *M*-* a< £2 >F(k l7 k,) 

(102) 



F(k i5 k,.) = ^2^ 2 '(q;y) e - <q - (x - x,) q 2 

q 

+[^+^7 S A 2'i' (q; y)(-*q)+27 E (A 2 2 ' (q; y')(*q) + ^'i(-q; y)(«0) e" iq - (x " 

q q 

[* k J + ^7 S ^ 21 ( q; y )(- zc i)+^7 S (A 2 2 ' (q; y')(-*q) + (q; y)(-*q)) e ~ 4q ' (> 



q q 



(103) 

q q q 

E F22 ( q; ") + ^7 S Asi fe^H^E A 2' 1' (q; *. j) 

q q q 

-7^7 E I A 22' (q; *. i) + ^2' 1 (-q; h 3) + -V (q; h 3) + (q; *» e- iq -< x - x ) 

(104) 

Fi 2 (q;y) = P 4 (q)P J (q)[-z(k t .q)q 2 A 2 ] (105) 
F 22 (q;y) =^(0)^(0)^] (106) 
F n (q;y) = ^(q^qH^.qXk^A 2 ] (107) 
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P 2 i(q;*j) = Pi(q)P J -(q)[i(k J -.q)q^A J ] 



^(q)= Hb + 



(k,.q) : 



Define, 



4m to 



A(q) = -k«.q+y 



^ 2 'i(q;^i) = 



2to 2 A 2 



A(-q)^(q) 
2m 2 A 2 



A 2 (q) 

. _ 2m 2 A 2 

^2i'(q;*>.?) - D 2 ( q ) 



2m 2 A 2 
A(q)^(q) 



2 A 2 



2F 



1 



A(q) A(-q) 



(108) 
(109) 

(110) 
(111) 

(112) 
(113) 
(114) 
(115) 



Define, 



L < E 2 >= t < E 2 > 
2 2 



2 A 2 



2F 



1 



2 to 2 A 2 



A(q) A(-q)y 2v ^VA(q) A(-q) 



1 



- <P 2 >= — 



V 



2 ™2 



(116) 



A 2 



?Ui(q) A(q)J + ^ I^VA(-q) + A(q)j e 



-iq.(x— x ) 



The above equation leads to a finite result as we may see below. 



(117) 



1 ,„ to 2 A 2 
-<P 2 > = 



V 



A(-q) AW/ VA(-q) A(q) 



1 



1 



V 



TO 



A 2 



— 2 



2V 



1 



1 



1 



A(-q) A(q) A(q) A(-q) 
1 1 



^A(-q) A(q) A(q) A(-q) 



sm[q.(x-x)] (118) 
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Similarly, wc may write for F(kj,kj) as follows, 

F(k k ) - — V 2m 2 A 2 q 2 iq , (x _ x ' } 

.., 1 ^ 2m 2 A 2 ... 1 v^/ 2m 2 A 2 .. . 2m 2 A 2 .. \ _ ia fx _ x ' h 

... 1 >^ 2m 2 A 2 ... 1 x^f 2m 2 A 2 . . . 2m 2 A 2 . . A _ ia , X _ X V 

(119) 

Now we would like to solve for F(kj,kj) and | < E 2 > in one dimension and 
three dimensions where the integrals are likely to be simple. 

B.l : One Dimension : 



, 2 to 2 A 2 /l 1\ 2 , , m 2 A 2 , J\ 1 . 



^2^2 ^ — 2ikjX ^2ikiX ^ ^2ikix 

2 kikj kj(^kj ~\~ kj^j k^i^kj ~\~ kj^ ~ k%k-i k^k^ 

jjj2 ^2 ^ ^ — 2ikj x ^2ikiX ^ ^ — 2ikj x 
' — I - ^ — \~ ~\~ }sf/72^X^] 

2 k^kj kj(^kj ~\~ kj^ k^i^kj ~\~ k^j ^j^j ^j^j 

(121) 

f(ki,kj) = {JijiU^-JjiiUau)) = - T dx e'to-^e^^Ffakj) 

*-* J-oo 

(122) 

The final formula for the a.c. conductivity involves evaluating the following 
integral. For muj « k F we have, 

J / 1—kF+Oj/vF rk F \ 

i?e[cr(w)] / dfcj + / dfcj [/(fc J (l+mw/4),i: J )+/(-fc 3 (l+mw/4): fc 3)] 

w \J-k F Jk F -u>/v F J 

~ — [/(Ml + muj/k 2 F ),k F ) + /(-Ml + mw/fcf.), M 

+/(-fcjr(l + muj/k 2 F ) 7 -k F ) + f(k F (l + mu/k F ), -k F )] (123) 

We would like to systematically evaluate this and show that in the zero frequency 
limit, the a.c. conductivity is proprotional to ui and vanishes for w = 0. The 
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uj — > limit is quite subtle. In particular, we may not set u> = at the outset. 
For then, f(k F , k F ) = oo. 

1 I" 00 mu, i "* 3 ^ 2 " J gn(x) 

f(k F (l+moj/k 2 F ),k F )^-- dxe^~ x e *f F(k F ,k F ) (124) 

^ J-oo 

F(k F ,k F ) = -i r ^-^{e~ 2ikFX - e 2lkFX ) sgn(x) 
2k F 



m 2 A 2 g^ikpx e ~2ik F x 



+ [ik F + ^—{-—^ + —^}sgn{x)Y (125) 



and, 



f(-k F (l + muj/k F ), k F ) « f(-k F , k F ) 

2m 2 A 2 l„l .- 2m 2 A 2 



- y / e 2lfc ^ e e F(-k F ,k F ) (126) 



F(-k F ,k F ) = - (2m 2 A 2 ) e - 2 ^* 
m 2 A 2 



+ — £3 — (— 1 + e- 2ifc ^) Sff n(a;)] 2 (127) 



1 Z" 00 : m 3 £\ 2 ^ sqn ( x ) 

f(-k F (l + mu>/k 2 F ),-k F ) w -- / da; e'F 1 e 'T F(-k F ,-k F ) 

L J -co 



(128) 

F(-k F ,-k F ) = F(k F ,k F ) (129) 



and, 

/(fc F (l + muj/k 2 F ),-k F ) rts f(k F , -k F ) 



- / d^e-^e V e^T F{k F ,-k F ) (130) 
F{k F , —k F ) = - (2m 2 A 2 ) e 2 ^* a: 
-[-*A!f + + e2ifcF3; ) *W*)] 2 ( 131 ) 
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f(k F (l + muj/k 2 F ), k F ) + f(-k F (l + mco/k F ),-k F ) 



(A 4 fc F m 5 ^ + (84 + A 4 m 4 )(164 - m 2 w 2 )) sm [A 2 m 3 u;/k F } 



4fc 4 f^(16fc 4 -m 2 cj 2 ) 
f(k F {l + muj/k F ), k F ) + f(-k F (l + muj/kp), -k F ) 



8A 2 k F m 2 



A 6 m 6 A 6 m 8 (61fc| + 8A 4 m 4 )w 2 



f(-k F (l+mu>/k%),k F ) = - 



+ 



96kfL 



+ 0[u 



(132) 



(133) 



fc 4 [A 4 m 4 cos[2A 2 m 2 /fc3,] + t (k F + iA 2 k 3 F m 2 + A 4 m 4 )sin[2A 2 m 2 /k F ]] 

~T A 2 m 2 (fc| + A 4 m 4 ) 

(134) 



These have been evaluated using Mathematical M . Immediately we see several 
problems. f(ki, kj) > for all arguments. Yet the above equations tell us that 
they are negative sometimes and sometimes even complex ! 

This means that something is seriously wrong with our formalism. Perhaps 
replacing V6*i • VAj by the average (which is zero !) was not such a good idea 
after all. 

What is as bad is that the final formula for the a.c. conductivity has a non- 
vanishing constant part (at least it is real !) this means that d.c. conductivity 
is not zero. This is clearly wrong. Maybe some readers of this preprint will offer 
to collaborate with the authors to fix this difficulty. 



B.2 : Two and Three Dimensions : 



In two and three dimensions, the integrals are substantially more complicated 
caused by the complicated angular parts. Work is in progress in collaboration 
with Shri. Chandradew Sharma it will be reported soon. Useful conversations 
with Prof. N.D. Haridass of IMSc. is gratefully acknowledged. 



6 Appendix C 

We would like to ascertain whether or not the zero frequency limit of the a.c. 
conductivity is the d.c. conductivity. First let us define d.c. conductivity. 
Consider a system of electrons coupled first to a uniform d.c. electric field. The 
interaction part of the hamiltonian may be written as, 

Hi = \e\ f d d x ^t( x ) ( Eext . x) ^(x) (135) 
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The expectation value of the total momentum of the electrons may be written 
in the interaction representation as, 



T S P(t) 

(P(*)> = d36) 

s = e -ij; %0 dt (137) 

The d.c. conductivity is then simply given by, 



Since, 



= ^1^< P <'»I < 138 » 



5E 



ex l 



r—i/3 , f , , , , , , 

S = -i J dt\e\ J d d x'^(x',t) x' ^(x',t) 



(139) 



X(t') = y AV^x',/) x ip(x,i) (140) 
<7d.c. = -* -^7 / * dt (TX(t)-P(tj) (141) 



Define, 



From the main text we see that the zero frequency limit of the a.c. conductivity 
is written as, 

a a . c .(0) = / dt a(t-t) 
Jo 

J dt J dt\ 0(t — ti) {T P(ti) • P(i)/ (142) 



? e 



m 2y 

Here P = k c^Ck- At first sight it appears that <Td. c . of Eq.( 141) is not 
equal to er a . c .(0) of Eq.( 142). For the two expressions to be the same we must 
have, 

X(t) = — / dhOti -ti)P(ti) (143) 



If we take the derivative with respect to t we find, 

!*(*) = (144) 

This is nothing but the definition of the momentum operator. It is a trivial 
kinematical result. Hence we may conclude that the zero frequency limit of the 
a.c. conductivity is in fact the d.c. conductivity. 
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